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We analytically perform the transverse momentum integrations in the real cor- 
rections to the longitudinal 72 impact factor. The resulting integrals are Feynman 



■ parameter integrals, and we provide a MATHEMATICA file which contains the inte- 

' grands. The remaining integrals are carried out numerically. We perform a numerical 

test, and we compute those parts of the impact factor which depend upon the en- 
ergy scale sq: they are found to be negative and, with decreasing values of sq, their 
absolute value increases. 



X ■ I- INTRODUCTION 

The NLO corrections to the 7* impact factor are calculated from the photon-Reggeon ver- 
tices for qq and qqg production, respectively. NLO corrections to the qq intermediate state 
involve the production vertex at one-loop level, rlV_^„:j, These virtual corrections have been 
calculated in |l[ As to the real corrections, the squared vertex |r^1'*_^gggp is needed at 
tree level; it has been computed in J3| and M for longitudinal and transverse photon po- 
larisation, respectively (cf. also [3, In Q we have combined the infrared divergences 
of the virtual and of the real parts, and we have demonstrated their cancellation. What 
remains to complete the NLO calculation of the photon impact factor are the integrations 
over the qq and qqg phase space, respectively. A slightly different approach of calculating 

n, 

the NLO corrections of the photon impact factor has been proposed in [Qj. Recently, the 



NLO calculation of another impact factor has been completed \^, the impact factor for the 
transition of a virtual photon to a light vector meson. 
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In this paper we perform, for the case of the longitudinal photon polarisation, the phase 
space integration in the real corrections. Our aim is to have, as long as possible, analytic 
expressions which allow for further theoretical investigations. The main obstacle is the (in- 
finite) integration over transverse momenta: in order to be able to perform the integration 
analytically, we introduce Feynman parameters. This will allow for further theoretical in- 
vestigations of the photon impact factor. In particular, the Mellin transform of the real 
corrections w.r.t the Reggeon momentum can be calculated. This representation (together 
with an analogous representation of the virtual corrections) will also allow to study the im- 
pact factor in the collinear limit and to compare with known NLO results; it can also be a 
starting point for the resummation of the next-to-leading logs(l/x) in the quark anomalous 
dimensions. 

Starting from the results of j^, we have to integrate a sum of expressions, corresponding 
to products of Feynman diagrams ^ that differ in their denominator structure. In order to 
introduce Feynman parameters we therefore split this sum and treat each Feynman diagram 
(or small groups of them) independently. This gives rise to divergences, which in the sum 
of all diagrams will cancel but show up in individual diagrams. The main task in the 
program of performing the integration analytically is the regular isat ion of these divergences 
in each individual diagram: for this we use the subtraction method. After carrying out 
the integration over transverse momenta we arrive at analytic expressions for each diagram. 
These expressions are convergent integrals over the Feynman parameters and the momentum 
fractions of the quark and the gluon. Since the expressions are somewhat lengthy, we do 
not list them in this paper but provide a MATHEMATICA code which we describe in 
the appendix. As a first application of our results, we carry out the remaining parameter 
integrals numerically and perform first tests of our calculation. In particular, we study the 
dependence of the photon impact factor on the energy scale which is entirely contained in 
the real corrections: it turns out to be in agreement with the expectations. 

Our paper is organized as follows. In the next section we will recall the structure of the 
real corrections and specify the expressions we have to integrate, following very much the 
strategy described in Q]. In the following section we introduce Feynman parameters and 
perform analytically the integration over transverse momenta. The main emphasis lies on 
the regularisation of the divergences, which appear in the individual diagrams. Finally, the 
numerical results are discussed. An appendix describes the files that provide the Feynman 
parameter representations of the individual diagrams. 

^ For simplicity we will in the following simply use 'Feynman diagram' rather than 'product of Feynman 
diagrams'. 
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II. THE REAL CORRECTIONS 

The starting point for the real corrections is the squared vertex |r^'i''_^^ggp which, starting 
from the process 7*+g QQg+Q for the longitudinal photon polarisation, has been calculated 
in jsf. The notations are shown in figO The momenta are expressed in terms of Sudakov 
variables: k = aq' + Pp + k±, i = aeq' + Pep + i±, with q' = q + xp, = —q^ and k"^ = —kj_. 
FigEl shows the Feynman diagrams. The product of two diagrams, summed over helicities 

q — k — £ q 




FIG. 1: Kinematics of the process 7* + ^ qqg + q' ■ 



Ai A2 A3 A4 




FIG. 2: The dia grams contributing to r^*^^^. 



and colors of the external particles, is labelled following the notation of the diagrams (e.g.: 
AB12 = A1B2). The expressions for the products AAij, ABij, ... that we use here are taken 
from |2j where, with one exception, we have used the same notation: in the present paper 
we combine, for the sake of simplification, the diagrams 4 and 5 in the product with diagram 
4, i.e. we set: 

AA44 + AA45 + AA54 — > AA44, 

^^5,^^54 ^ 0, 
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ABu + AB45 + AB54 — > ABi4, 
AB45,AB5A — > 0, 

and similarly for BA and BB. All products AA etc. have (transverse) dimension —6; 
in order to deal with dimensionless expressions we multiply with (Q^)^. Let us introduce 
dimensionless variables: 

k, £, r, A, ^ y^(fc, £, r, A, (1) 

(the meaning of A and Sq will become clear in a moment). From now on we will use only 
these new variables. We also define the abbreviations 

ai = a , 

ai = {1 — a — ai) , 
^2 = (1 — a) , 
0:2 = (a + ai). 

Furthermore, we introduce the label JB, in order to denote a generic product of amplitudes 
iAA,AB,BA or BB): 

A3 = ^-^(^^11 or AA12 or AA21 ... or ABu. . . .). (2) 

It is then easy to see that these jB^s are dimensionless and only depend on the dimensionless 
momenta k, £ and r. In Q we have included a part of the qqg phase space measure in the 
definition of the ^'s, in order to simplify the expressions below. Finally, |r^°''_^gggp will be 
proportional to the sum of all the 



The procedme of arriving at finite NLO corrections to the 7* impact factor has been de- 
scribed in 



Let us briefly review the main steps. |r^?_^^qgp has to be integrated over 
the qqg phase space. Before doing this integration, two restrictions have to be observed. 
First, we have to exclude that region of phase space where the gluon is separated in rapid- 
ity from the qq pair (central region); this configuration belongs to the LLA and has to be 
subtracted. To divide the qqg phase space an energy cutoff sq is introduced. This energy 
cutoff plays the role of the energy scale: when combining the NLO impact factor with the 
NLO BFKL Green function it will be important to use the same scale in all pieces. Since 
the virtual corrections to the NLO impact factor are independent of Sq, all dependence on 
the energy scale resides in the real corrections. As a result, the calculations described in this 
paper can already be used to study the sq dependence of the NLO impact factor. Next, we 
need to take care of the infrared infinities. The divergences in |r^°''_^g^gp due to the gluon 
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being either soft or collinear to either of the fermions are regularised by subtracting the 
approximation of the squared vertex in the corresponding hmit. These expressions are then 
re-added and integrated in D = 4 — 2e space-time dimensions, giving rise to poles in e, which 
drop out in combination with the virtual corrections and to finite pieces. The subtraction 
of the collinear limit requires the introduction of a momentum cutoff parameter, A. The 
final NLO corrections must be independent of this auxiliary parameter. In our numerical 
analysis we will perform this important test. 

According to the full NLO corrections to the 7* impact factor have the following form: 

finite 2$^°'' f 



_ ^(l,virtual) 



"7 



■7' 



(47r)2 {"^'^ ^ 



l3o\n— + CFHr' 



J dc( J ^^^y fc) |Ca [In^ a{l - a)so - In^ M^] 



+ 2C 



F 



8 - 3 In a(l - a) A^ + In^ + In^ 



a 



1 — a 



(l,rea/) 



finite 



Ca ^ 



{l,real) 



finite 



Cf 



(3) 



The terms multiplying stem from the UV renormalisation which belongs to the virtual 
corrections and does not need to be discussed in the present context. The finite pieces which 
are left after combining the IR singular pieces of the virtual and real corrections above are 
given in the second and third lines of eq.Q. I2 is proportional to the squared LO photon 
wave function: 

with D{k) = fc^ + a(l — a). Finally, the squared mass of the quark-antiquark pair, M^, is 
a function of a, fc, r, and it is not important for our present analysis. 

Let us now focus on the last line of eq.Q. In terms of the products of single diagrams, JB, 
these finite contributions to the real corrections yi read: 



^ ^il,real)f-^'^ f^VW^ f ,2-2e , ,2-2e , f' [''''da, 




?soft - (^q - ^q,.o/t) e(a,A - 

- (^q - JS^,soft) e(a,A - \t'\) i, (5) 
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X 



Ca 



J3 - JSsoft - (^con - ^cen,soft) 0(|^| - "^V^) 



- r — V^cene(|£|-a,v^)l (6) 

The evaluation of these two expressions is the issue of this paper. The sums extend over 
the Cp and Ca parts of the ^'s, respectively. Note that the sum over all JB^s is essentially 
|r^?_^gqgp. Whereas the full sum is finite, the individual contributions ^'s, without further 
modifications, would be divergent: our task, therefore is to render them finite, diagram by 
diagram. This is the content of eqs.© and In order to indicate the different limits 
where the integral of JB diverges we use the following subscripts: 

soft := ae ~ |€| ^ 

q := \e'\ 

q := \£"\ 
cen := — 

^X,S0ft:=[^x]3off (7) 

Here £' i" denote the two coUinear vectors 

£' = £+^k , £" = £-^(k + r). (8) 

The parameter A defines a cone around the coUinear directions, specifying thereby the region 
of the subtraction. The exclusion of the central configuration (cen) is realised by subtraction: 
only the rapidity region a£ > \£\/y/so is counted as a contribution to the impact factor. 



III. THE INTEGRATION 



Our aim is, in eqs. ()5l6|) . the analytic integration over k and £. To this end we introduce 
Feynman parameters. Since different Feynman diagrams provide different denominators, 
we have to deal the JB^s independently rather than in the sum of all (as an alternative, 
an attempt the find common denominators would lead to expressions that have too lengthy 
numerators). As indicated in eqs. ()5l6p . in each single diagram JB we have to combine the full 
expression with subtractions which are obtained from certain approximations (soft, cen,...), 
given by eq.©. As a first step, we analyse the divergences of all diagrams. The table HI shows. 
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TABLE I: The divergences of the diagrams AAij and ABij 



which diagram diverges in which hmit. Both in the left hand and in the right hand part of 
the table, it is the second column which lists the divergent limits which in eqs. (|5l6p require a 
suitable subtraction. The third columns contain an additional divergence that appears due 
to the separate treatment of the diagrams; we shall take care of them by making a further 
appropriate subtraction form each j^. Let us go through these subtractions in some detail. 
Each product JB has the general structure 

^ = —-^-— , D^ = D^{k,£,r). (9) 

The denominators D„ are given by (here we use definitions which slightly differ from those 
of B): 



Di = at [aiai + ai{k + i + rf + ai{k + rf] + aiotif (10) 

^2 = aias + (fc + r)2 (11) 

D^ = {ai£-at{k + r)f =a\i'^ (12) 

Di = {ail + at{k + t)f =al(!'^ (13) 

= an (aiOi + ai{k + if + ai{k + rf) + aiai{i - rf (14) 

Dq = ae [aiai + aik"^ + ai{k + if] + aiai i? (15) 

^7 = aia2 + (fc + £)' . (16) 
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Note that, in some of the ^'s, two of the four D's in eq.® may coincide. The numerator 
Z is a polynomial in scalar products of the (dimensionless) transverse momenta: 

Z = X^k"^ + X2£'^ + X3k£ + X^kr + X^ir + Xqv'^ + Xt , Xi = Xi{a,ae). (17) 

(without our choice of dimensionless momenta, (^, the last term, A7, would have been 
proportional to Q"^). Let us now consider one single JB. Using 

0.^1, '5(1 -E. A) 



we introduce the Feynman parameter representation: 

J3 = Q 5(1 - 5^ A) f, (18) 

The j3i carry the same labels as the denominators Di they belong to, and the sum in the 
denominator of (fTHj) is understood to extend over those indices that occur in the diagram. 
//3 stands for products and powers of the I3's in the numerator which appear in case of some 
of the four D's being equal. For example, the product ABi-^ reads 

Z , (5(l-EA)/3i^ ..Q^ 

Di V2 V3 Jq [PiDi + /32-D2 + P3-D3J 

In order to determine the Feynman representation of the approximations of JB, we first 
consider the numerator and the denominator, in eq.Q, in the appropriate approximation; 
we denote the result by a subscript. Then, the Feynman representation is introduced as just 
described: 

J3, = 6[ r[rfA,5(l-EA)//3 rv^/hi x = g, g, cen, ... (20) 

Each limit in the original momentum space representation (in the {k,£,ae} space) un- 
ambiguously corresponds to a 'corner' in the {k,£,ae, f3i} space of the Feynman parameter 
representation. By our way of calculating the Feynman representation of ^soft? ■ ■ 

((201), ^6 ensure that the cancellation between JB and its approximations, which originally 
was formulated in the momentum space Feynman amplitudes, remains valid also after the 
introduction of Feynman parameters. The prescription (pUj) in particular ensures that our 
original expressions for the real corrections, eqs.(jni EI), are exactly translated into the Feyn- 
man parameter space. 

We now perform the integrations over k and £. According to the different integration regions 
that occur in eqs. ()5f(ij) we are faced with three different types of integrals. Let us discuss 
them in some detail. For convenience, we generalize to 2 — 2e transverse dimensions. 
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Let us consider the integration J (P~'^^k (P^'^'^l JB. The ^-integration extends over the whole 
space. After the introduction of Feynman parameters, eq. ()18|l . we have to integrate 

Ji := I £-'^k d'-''£ (21) 

The generic structure of transverse momenta is: 

J = f £-2ej^^ ^2-2e^ ^1 + ^2 + A3 k£ + A4 fcr + A5 Iv + Ag + A7 .^2) 

^ ' [rk'^ + p£^ + 2-fk£ + 2r]kr + 2ij£r + Dr'^ 

The coefficients Aj in the numerator depend upon a, ag. The coefficients in the denominator, 

r, p, 7, are functions of a, ae and of the Feynman parameters f3i. Dr depends, in addition, 

on (and on Q^, which has been normahzed to 1). In order to ehminate, in the denominator, 

the mixed scalar products, we perform the shifts: 

7 1 
k = k' --£' - -^^{PV - 



J, = / s-^-k- <p---e ,-^5^ + ^.5' (24) 




We arrive at: 

.h= [ ^2-2. fc/ ^2-2.;,/ ^ 

where we have introduced the definitions 

u = - detfi , R = Dr- -r^iTfp + r/i^ - r//i7)r2 , Vt = \ ' ' 1 . (25) 
r detiz 

In the numerator of ()24j] we have already dropped the mixed scalar products {k' k' r', £' r') : 
when performing, in (j221); the shifts, such scalar products appear. However, after angular 
integration their contribution vanishes. Note that except for A = Xi all coefficient func- 
tions appearing in eq. lj^ . in general, depend on the Feynman parameters Pi and on a,ae. 
Performing the momentum integrations we get: 

j, = [ d'-'^k d'-'^£ + ^5 

J [rfc2 + u£^ + RY 

e)r(l-2e) - + - +r(2 + 2e)- ^. (26) 



Note that the integral converges for e = 0. We will need the general result (e 7^ 0) later. As 
our final result for the Feynman parameter representation of JB we define a quantity X by 

j d^kdHj3 = 6j^' l[dp,6{l-J:P^) U Jil=o 

= J n^A5(i-EA) X. (27) 
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For each diagram J6, the function X is given as MATHEMATICA code, described in the 
appendix. 

In order to integrate the soft (£ ~ ^ 0) approximation of JB we consider the numerator 
and the denominator of JB in the soft approximation and introduce the Feynman parameter 
representation according to ((201) • We then carry out the momentum integration exactly as 
just described ending up with an expression which we call Xgon'- 



[ d'kdHJ3,oft = Q [ n^A5(i-EA) [d'kd'i 



soft 



* I softJ 



(28) 



It is instructive to see how the soft limit in the momentum space translates into the Feynman 
parameter space. The corresponding region in the Feynman parameters space is determined 
by the behaviour of the denominators Di in the soft limit. One obtains EA-^j|soft from 
EAA by 'weighting' {/32,/?5,A} ^ {A, A, /?7}p^ {A,^,-^,"^} ^ {Pi, f3e,£,ae}p, and 
then expanding around p = and by keeping only the most divergent term. The prescription 
for the /5's remains unchanged after the integration over £ and k; the soft region is therefore 
specified by: 



soft : 



Hence, we can obtain Xgoft either in the way given in eq. 
in the limit (indicated by the subscript): 

d'^k d'^i JSsoit 



« A, A. (29) 
or by considering X feq. ()27p ) 



1 d^k d^e 









A, A, A < A, A, a^, £ < A, A 



nc?A'5(i-EA) [X] 



(30) 



Note that we must not apply the approximation ()29|) to the argument (5-function in eq. ()30|) 
(otherwise, the translation of the real corrections from the momentum space to the Feynman 
parameter space would be incorrect). This has the following consequence. After choosing 
any of the /9-integrations to be done by means of the 5-function we have 



j d(3k n d(3, 6{l - E A) {X - Xsoft) = j\[d(3, (x- X,oft 



(31) 



with the tilde symbol on the rhs indicating that (3k is expressed in terms of the other j3i. 
Since the argument of the 5-function is not approximated, Xsoft is not equal to X taken in 
the limit (|^. However, both Xsoft and X become equal in this limit. 
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• JBcow 

Next we turn to the coUinear approximation of jB. We use the label 'coll' to denote the 
two collinear limits which are defined as either \£"\ or \£'\ (cf. eqs.Q and 
dHl))- Following the notation introduced in eq. (j20|) we use the collinear approximation of the 
numerator of JB and of the D's. The region of integration is restricted to a cone around the 
collinear direction: 



J, := / d'-'-k d'-'-e e{aeA - |4|) ,4 > = i\ t. (32) 



«lcollJ 



Starting from eq. (j32j) we first express i in eq . (j32j) through if.. The list of the D' s in eqs. fjlUp 
- (fTBj) shows that in the limit (! 0, in leading order, all Di become independent of 
except for D/^ which is proportional to The same holds for the the second collinear limit, 
with Dj, instead of D^. The denominator in eq. ljH^ therefore depends on only via £^ and 
the integrand is parametrised as 

J2 = y d^-'^'k rf2-2^4 0(a^A - 141) 

^ Ai k^ + X2£l + A3 fc4 + A4 fcr + As IcT + As + A7 
^ [Tk^ + p£l + 2r]kr + Drf 

Note that, since we are using i^ instead of £, the coefficients in eq.(|33|) are not the same as in 
eq.()22|). taken in a collinear approximation. It is only for keeping the notation as simple as 
possible that we do not introduce new names for the coefficients. In order to diagonalise the 
denominator we only have to get rid of kr; the momentum £c, therefore, does not participate 
in the shift of momenta (j23|) . and the limits of the £c integration remain the same. Applying, 
to (jnSl), the shift with 7 = /i = 0, we arrive at 



TT^A^ r 1 1 A2 1 2i? + A2^ 



6tujR{R + A^) 



with : A^ = uj{aeAy. 



(34) 



Below we will show that only the result for e = is needed. In analogy to eq. ()27|1 we define 
Xcou to be 6//3J2 after performing the integration over the transverse momenta: 

j d'k d'£^,,n e(«,A - 141) = 6^' l[dp, 5(1 - E A) //3 J2l=o 

= J n^A5(l-EA) ^cou. (35) 
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In the collinear limits all Dj, except (D^), become independent of £c- In the f3i space 
these limits therefore read: 

collg: (3i^4 < Pa (36) 

collg: < Ps. (37) 

However, in momentum space the integrals of .^coii and have different limits of integra- 
tion. It turns out that, as a consequence, Xcou, in the limit (jHU)) or (jSTj), can be written as 
a sum of two parts. The first one coincides with the collinear limit (either (jHUj) or ()37|)) of 
X. The other term depends on A, and vanishes in the collinear limit ()36|) or (j?fj) . 

Next we need to consider the soft limit of the collinear limit, the collinear-soft approximation 
"^coiLsoft- In momentum space, this limit is calculated by taking, in JBco\\{k,£c,ae), the 
additional limit £c ~ ^ (cf. the definition of £c in It is then integrated exactly in 
the same way as just described for Xcou, 

J (fk dH ^coii,soft e(a,A - 141) = j riA '5(1 - E A) Xcoii,soft. (38) 

Turning to the collinear-soft limit in the (3i space, one finds that the soft limits of the Di 
and of the -Dilcou coincide. However, the ^c-integration of ^coii,soft in (|HH|). in contrast to 
the integral of ^soft, has an upper limit of integration. But one can show that due to the 
ardependence of the integration limit a^A in eg. ()38j) the position of the soft limit in the 
subspace (of the {^c^^f. A} space) is unchanged by the 4 integration. Also in case 
of the integration of ^soft (without upper limit on £) the soft region is unchanged by the £ 
integration, see eq. (j3(J|) . The soft limits of X and Xcou, therefore, are located in the same 
region of the jSi space: 



softofcoU: A, A, A < /?i,/?6,a£ < A, A- (39) 



The third type of integral appearing in the real corrections, eq.® and ©, deals with the 
central approximation JBcen^ defined by the limit — > 0: 



J, := I d'-'^k d'''^£ Q{\£\ - aeso) (40) 



icenJ 



Its general form reads: 

J3 = /" d^-^'k d^-^'£ ei\£\ - a^so) 



Ai fc2 + A2 ^2 + A3 fc£ + A4 fcr + A5 €r + Ae + A7 , , , , 

X J . (41) 

[rk'^ + p£'^ + 2-fk£ + 2r]kr + 2i2£r + Dr] 
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In order to keep the region of the ^-integration as simple as possible, we do not perform 
any shift in this will leave us with angular dependent terms of the form £r, both in the 
numerator and in the denominator. We only perform the shift eq. lf^Hj) of k, given by eq. ()23p . 
with £' being expressed through £ and r: 

1 



(42) 



and we obtain 



J3 = y d^k Q{\£\ - a^v^) 



Ak"^ + Be^ + C + E£r 
[rfc2 + uj{£ + Ar)2 + R]^ 



[-Ca + -Cb + -Cc + -a r Ce 

otujK \ t uj K K 



(43) 



The coefficients u and R are given by eg . (1251) . expressed in terms of the coefficients in the 
denominator of eq. ()41|l . Again, we only need the case e = 0. The coefficients in the second 
line on the rhs of (jlH]) are given by 



Ca 
Cb 



2 1 K 



1 \\^ + R 
2 



R 



so-{R + \^ 

So {R + A^) - 3i?2 - 



1 

Rk 



Cc 
Ce 



k3 

2 1 K 



(R - X'){R + Xy + So {R' + X- 6RX' 



2 1 Kl 



so + {R- A^) 
so + {R- A' 



2R 



2R 



So {R-X^) + {R + XY 
So {3R-X^) + iR + X'f 



where we have defined : 



(44) 



So = uJ[ae^JSo) , 



2„2 



and: k = {so + R + X^ - Asq X"^. 
In complete analogy to the previous integrations we finally define 

J d^k dH ^ccn e(|£| - «,v^) = j riA 8{i - E A) 



(45) 



Turning to the soft limit of the central-soft approximation, we again start in momentum 
space, introduce Feynman parameters, perform the momentum integration, and arrive at 
Xcen.soft, in analogy to eq.ip^. Due to the different limits of integration of JB and JBcem 
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only a part of X^en matches the central approximation of X, the rest depends on Sq and 
vanishes in the central limit. Similar to the collinear approximation, the ar dependence of 
the lower limit on the I integration has the effect that the momentum integration does 
not change the position of the soft region in the {a^, [3i} space. The soft limit is located in 
the region: 

soft of cen: ^2,^5,^7 < < /^s, /?4, /^e- (46) 



As we have mentioned before, when treating individual diagrams, JB, additional logarithmic 
divergences appear which do not show up in the sum. They occur in the following limits: 

e-uv : 1^1 ^ 00 

is : \£\ ~ y/a^ — » 

\t : £^r. (47) 

Tab U contains, in the third columns, a list of those diagrams where these divergences appear. 
These divergences lead to additional subtractions which we have to describe in some detail. 
Let us start with the Cp parts. They do not contribute to the central limit. The only 
additional divergence is of the type i — nv; in this limit, all denominators Di |TO|l - |T6| 
are proportional to i"^, except for D2 which is independent of £. Any ^ containing in 
the denominator and a term proportional to £^ in the numerator will therefore, in the limit 
£ 00, be proportional to l/£'^, leading to an UV divergent £ integration. Subtracting from 
such diagram its approximation in this UV-limit, JBi-uv, would cancel the £—nv divergence. 
However, the £ integration then becomes IR divergent instead, since JB^-av ~ l/£^. We 
therefore define our subtraction term in the Feynman parameter space where it is easy to 
avoid this additional IR divergence. Let us demonstrate this at a simplified expression 
A = l/[(c + £'^)D2]. c is a constant and D2 does not depend on £. We introduce the 
Feynman representation and integrate over £ (assuming already a term that cancels the 
£ — uv divergence): 

' d£' A = r de- \—- = [ d(3,d(32 r d£' 



c + £^)D2 J "Vo m£^ + c) + f32D2[ 



2 



The divergence at £ — 00 appears in the limit f3i ^ after the integration. The natural 
subtraction in the Feynman parameter space therefore reads 







dPi dp 



(3,{(3,{c-D2) + D2) (3,D2 



(49) 
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As an alternative subtraction term we could use the £ — nv limit of A, which is 1/[£^D2], 
in Feynman parameter representation. It is obtained by taking the integrand in the second 
line of eq. ()48|) in the limit jSi —>■ excluding the argument of the (5-function from this 
approximation. Carrying out the (32 integration this results in 

1 

(3iil-(3i)D2 

The additional divergence of this term at /5i ^ 1 corresponds to the limit £ ^ in the 
momentum space and can be avoided by using the subtraction eq. ljl^ . 

We therefore regularise the £ — nv divergence of our diagrams in the following way. Due to 
the behaviour of the Dj's in the £ — uy limit mentioned above, the region £ —>■ oo translates 
to the region /5j^2 in the Feynman parameter space. Each £—uv divergent diagram has 
a D2 in the denominator. After the introduction of the Feynman representation and after 
carrying out the momentum integration we always perform the (32 integration by means of 
the (5-function in case of these diagrams. The subtraction term is then determined by taking 
the limit /3j^2 -C 1 in 2 — 2e transverse dimensions: 

j d?-^'k d^-^'£ J3= / rf/?2 n ^/^^ <5(i - /?2 - E A) ^ = /" n ^A X 

i i 

= j\[d(3, (x- + j\[d(3, (50) 

i i 

with: X^_uv = X 

To perform the momentum integration we needed to calculate the integral in eq. (j26p for 
e 7^ 0. The tilde in eq.fl^Uj) indicates that (32 in X is expressed in terms of the other (3i. We 
perform the (3i integration in the re-added piece analytically, obtaining an e-pole and finite 
terms. Summing up the re-added £ — uv subtractions from all diagrams we find that the 
poles cancel, as expected. 

We can now write down the final result (after momentum integration) for the Cp part of 
the real corrections (0). Using the definitions (jTfj) . (jHSj), and (jHHj) we obtain 



7' 



finite 

(l,rea() 



f da [ [ Yld(3J{l-EP^) 

^ Jo Jo «^ Jo 



Cp (27r)6 

— Xsoft — X^-uv — {Xq — Xg^soft) — {Xq — Xq^soft)^ 



Cp 



IT 



2 



+ ^{3n'-28)}. (51) 



16 



The superscript indicates that only the Cp part of the bracket is taken, excluding the color 
factor Cp- To simplify notations, we have written the delta function 6{1 — ^Iso for the 
X£_uv term. However, as discussed above, P2 in ^i-uv is understood to be already expressed 
in terms of the other j3i. We have set e = in (jKT|) since the jSi integrals are convergent by 
construction. This is why we needed to calculate the integrals over JBcoW and JBcen only for 
the case e = 0. The term in the third line of eq.()51|) results from the sum of the re-added 
£ — nv subtractions. 

Now we turn to the parts of the diagrams. We have to deal with collinear divergences. 
However, according to 'sj the collinear approximations of all Ca parts sum up to zero; we 
can therefore just subtract from the divergent JB's their collinear approximation without 
re-adding it (as in the Ca case, we restrict the integration to a cone since otherwise it would 
be UV divergent). 

According to tabUl in the Ca parts we encounter all three types of additional divergences. 
As to the is(intersoft)-divergence, one finds in momentum space: 

^ - J3,oft)is = 0. 

Ca 

We can therefore regularise this divergence by subtracting, from {JB — ^soft), its is- 
approximation, without re- adding it. In the appendix we provide the combinations 
{JB — ^soft)is in the Pi space, (X — Xsoft)is, which can be obtained from X — Xsoft by 
taking the is-limit (before making use of the ^-function): 

is: (32, (35, (37, ae < (3i, (3^, (3^, (3^. (52) 



In the Ir-limit {£ r), it is the central approximation of certain diagrams that diverges. 
Let us recall the structure of the AA and AB diagrams in the central limit — > (see jsf): 

af al Ca 1 



*J I cen 



AB 



(£-r)2 D^k 

— af al Ca 



D 



1 



U I cen 



(£-r)2 D{k + r)D{k 



where D{k) = axa^ + fc^. The matrix Dij reads : 



/ -i ^\ 






_ 1 

2 
r2 






_ 1 

2 

r! 

£2 



n _i ri 

^ 2 £2 




2 







0/ 



(53) 
(54) 
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The denominators in eqs.()53p and ()54|) correspond to : 

= a^a,{i - rf (55) 
D7U = D{k + £). 

Eqs. fj53|) and ()54j) show that, if a diagram contributes to the central hmit, its approximation 
in this hmit is divergent as £ — r. However, for a given pair (i,j), eg. ()53|1 and eg. ()54|1 become 
egual with opposite sign as € — > r, since -Dylcen ^ -^2lccn- regularise the uv divergence 
we therefore choose a common set of {Pi} for AAij and ^iSjj in the following way: 



The momentum integration is performed in the same way as for all other diagrams. For the 
Pi integration, we consider the sum {AAij + ABij). Our parameterisation ensures that 
the Ir-divergence (in the Pi space) cancels between the two diagrams. All what we have said 
about the Ir-divergence in the AAij and ABij in an analogous way also applies to the BE 
and BA diagrams. The parameterisation of BBij is given by eg. with D2 being replaced 
by D-j. Note that it is only those diagrams that are of type AA or BB and contribute to 
the central limit which need this special parameterisation. 

The matrix Dij reveals another potential divergence: ^^islccn -^-^zslcen obviously 
divergent in the soft limit. Their soft approximations which we will subtract are proportional 
to l/£^ giving rise to an UV divergent £ integration. However, one can show that this 
divergence is automatically cancelled by ^^15 1^^, and ^^25] soft' respectively. 

Following (jH)), we have to subtract the central and central-soft approximations in the region 
where \i\ > aiy/so- Egs. fl53|l and ()54|1 show that the central approximations of the diagrams 
14,24,41,42,44 are also £ — uv divergent. In order to regularise this £ — uv divergence in the 
same way as discussed above one would need to calculate an .^-integral with non-zero lower 
limit in arbitrary dimensions. This can be avoided since, according to the structure of Dij, 
24 41 42 44 -^cen = 0; wc are therefore free to subtract the central approximation for these 
diagrams in the whole momentum space {\£\ > 0) translating to: X — Xa^=Q. The £ — uy 
divergence then is regularised by subtracting the €— uv approximation of this combination, 
(X — Xq,^=o)£-uv, after doing the P2 integration by means of the 5-function. Since the central 
and £ — uv approximations commute, this is egual to (X^^uv — ^i-nv\ai=o)- We re-add these 
subtractions, carry out the remaining /3j-integrations and find the expected cancellation of 
the e poles. 
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Our final result for the Ca part, defined in eq.®, then reads: 



Unite 



(l.reai) 

1 rl-a J^. 



^ f da [ [ nrfA5(i-EA) 

1^ AA Jo Jo ^0 



12,14,24,44 



-'^ — -^£-uv — {X — X£_uv)a^=0 ~ -^soft — (-^coU " -^coU.soft) 

^ Ca 

X — Xsoft — {Xcen — -^cen,soft) — {X — Xsoh)is ~ {Xcoll — -^coU.soft) 



-E r — ^^rfA5(l-EA)x,^^^ 

^ Jo Jl-a Oii Jo 



Ca 

2 



7r2 



- (Stt^ - 28) ^ . (57) 

The summation in the second line is understood to extend over AA, AB, . . . and to include 
also 21,41,42; by 'rest' we mean all other Ca diagrams. The * on the sum indicates that 
in case of an Ir-divergent diagram we are not allowed to separate the sum of AA and AB 
[BB and BA) diagrams . The term in the last line again is a result of the re-added £ — uv 
subtractions. 

Let us summarise the procedure of deriving, from the A3's given in the FejTiman para- 
meter representation of the real corrections written in eqs. fl51|) and (jFfj) . The results of P|, 
repeated in our eqs.® and (jHI), define, for the sum of all diagrams, the necessary subtractions 
of coUinear and soft divergences and of the central region. We then turn to individual 
diagrams JB: we start from momentum space and calculate the various approximations 
required in (jSj) and ©• We then introduce the Feynman parameter representation, written 
in (jl8j) and (j2(J|) . retaining still the transverse momenta to be integrated over. In case of the 
collinear limit we change the integration variable £ to Ic- In the integrands, we explicitly 
write the squares and scalar products of transverse momenta. For each approximation we 
then obtain two sets of coefficient functions, {Aj} in the numerators and {r, p, 7, . . . ,Dr} 
in the denominators of the integrand, respectively. From the latter set we calculate the 
quantities uj,R,\ in and (jl21), and we define the appropriate shifts of the transverse 
momenta. Applying these shifts to the numerator gives the coefficients A, B, C, E. Now, we 
do the integration over transverse momenta, observing the necessary limits of integration. 
The results, eqs. fpB|) . (fH^. and are expressed in terms of the coefficient functions 
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obtained before. Using the definitions ^I7\i. and (jl^j) we are finally left with a set of 
expressions X, Xcen, -^cou, -^soft, • • • which are functions of the momentum fractions a,ae, of 
the Feynman parameters of the cutoff parameters A and Sq, and of the virtuality of 

the t-channel gluons, r^. After calculating X£_uv and (X — Xsoft)^^ we have got everything 
that in eqs. (jSH EZ|) is needed for a finite integral representation of the single diagram, JB. 

We have implemented this procedure into a MATHEMATICA program, and we have applied 
it to all JB's. The result are analytic expressions for all the X, Xcen, -^cou, • • • needed in 
eqs.(|^T|) and (fS7j). They are contained in a set of files described in the appendix. 



IV. NUMERICAL RESULTS 



For the final integrals in eqs. ()51l57p we have carried out the remaining integrations over 
a, a£, (3i numerically. We have done this for each JB separately, except for the Ir-divergent 
diagrams, where we have considered the combination AAij + ABij and BBij + BAij, re- 
spectively. We have used the Monte-Carlo routine VEGAS. For each diagram the integrand 
is given in a FORTRAN code which is described in the appendix. As a result, we have 
obtained values of 

Unite 



as a function of r , A and Sq. 



To be precise, we start from the expression Q of the NLO corrections to the 7* impact factor, 
and we combine the finite real corrections in the fourth line with the A and So-dependent 
pieces in the second and third line: 



with 



(1) 



^{l,virtual) 



+ j da 

'0 



finite 2$ 



(0) 



{Any 



/3oln— + Cf ln(r2 



{Any 



-31na(l-a) + ln^ + In 



a 



1 — a 



real 



(l,reai) 



= Ca 

+ Ca / da 
Jo 



finite 



Ca 



+ Cf <l>;t' 



(l,reai) 



finite 



(0) 



CF -^^-(#^^^ 



{Any 



■X2(Q;,fc)ln a{l — a)sQ 



e'ei {^A + $F + Aa + A,„^ 



(58) 



(59) 
(60) 
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In eq. fj60p we have introduced short hand notations for the four terms in eq.()59p. The LO 
impact factor is given by <l>^°2 = / d^k dal2{a, k), with I2 being defined in eq.(jl}. The 
integration of Aa, A^^ is straightforward and we obtain as a function of r^, A and Sq. 

In our numerical analysis we restrict ourselves to two questions. First we investigate the 
A dependence of the real corrections. The cutoff parameter A specifies the region of the 
subtraction of the coUinear approximations. Since these subtracted terms are re-added, the 
NLO corrections to the impact factor must be independent of A. This must also be true 
for ^^^f'', since it contains all A dependent terms. A^^ in eq. (p]|) does not depend on A, so 
+ Aa and must be A-independent separately. 

Second, we study the sq dependence of the real corrections. Since the entire sq dependence 
of the NLO corrections is contained in the real corrections and in the finite term fcf.eqs. (j59t 
IBUj) . the results obtained in the present paper already allow to study the scale dependence of 
the NLO 7* impact factor. Whereas the complete scattering amplitude, involving the NLO 
impact factors and the NLO BFKL Green's function, has to be invariant under changes 
of So, the impact factor alone is expected to vary. Since a decrease of Sq in the energy 
dependence {-^)'^ will enhance the scattering amplitude, the combined Sq dependence of the 
impact factors and of the BFKL Green's function has to compensate this growth. We find 
that the So-dependent part of the NLO corrections to the 7* impact factor has the sign 
opposite to the LO impact factor and, in absolute value, becomes more significant when sq 
becomes smaller. As a result, this part of the NLO corrections, in fact, tends to make the 
LO impact factor smaller. 

Fig El shows the dependence of the A dependent parts of $5^,1"' on the momentum cutoff A for 
= So = 1 and = 3. According to eq.(|59|) the only Ca term in ^1^*°-^ is ^a- It includes 
A dependent terms, since, as discussed above, we also performed collinear subtractions for 
the Ca parts of the diagrams. Fig|21 shows that, in fact, the A dependence is very weak. As 
to the Cp terms, $f turns out to be proportional to In A. This growth with A, however, 
is fully compensated by Aa. Note that $i? is a sum of many Feynman diagrams, whereas 
Aa ~ In A. The compensation of the A dependence, therefore, represents a rather stringent 
test of the calculation of the 7* impact factor. 

Next we address the dependence of the NLO 7* impact factor on the energy scale So. The 
full LO and NLO impact factor can be written as = g'^^^ + 9^^^^* ^ where g"^ = Anas- 
Since, at the moment, we only know the real corrections, we compute, as a part of the full 
answer: 
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FIG. 3: The dependence of the real corrections on A (with = sq = 1) 

We set e^Cj = 1. For the photon virtuahty we choose, as a typical value in 7*7* scattering, 
= 15 GeV^. This choice affects the strong coupling, a<(((5^) = 0.18 or (7 = 1.5, and, 
through (^, the normalization of and Sq. Fig0 compares $1^, to the LO impact factor 
g"^^^^} as a function of at different values of sq. In agreement with gauge invariance, 
both the LO impact factor and the real corrections vanish at = 0. The ratio of and 
g'^^^^} is shown in figEl The real corrections are negative and rather large. The overall 
magnitude is not of much significance, because we have considered only a part of the NLO 
corrections. More important is the fact that, in absolute terms, $1^, becomes more significant 
for smaller values of Sq. Since we included all Sq dependent terms in this implies that 
the 7* impact factor tends to become smaller with decreasing sq. This behaviour goes in 
the expected direction. 
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FIG. 4: , at different different values of sq 
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FIG. 5: The ratio <^>'^, / {g^ '^''^} ) at different sq 



V. CONCLUSION 

In two previous publications 0, Q| , the real corrections to the 7* impact factor have been 
computed, and in P| infrared finite combinations of real and virtual corrections have been 
obtained. The resulting integrals still contain the integration over the transverse momenta. 
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In the present paper we have carried out this integration for the real corrections, restricting 
ourselves to the case of longitudinal photon polarisation. To allow for further theoretical 
analysis we have introduced Feynman parameter representations, and we have performed 
the integration over the transverse momenta analytically. Our results are finite Feynman 
parameter integral expressions for each diagram. These expressions which may serve as a 
basis for future studies are presented in a computer code, which we describe in an appendix. 

In a first exploratory study we then have evaluated these integrals numerically. Two ques- 
tions have been addressed: first, we have shown that the NLO corrections to the impact 
factor are independent of the parameter A specifying the region of the coUinear subtraction. 
Second, we have studied the dependence upon the energy scale sq. A physical scattering 
amplitude (e.g. for the 7*7* scattering process), when consistently evaluated in NLO accu- 
racy, has to be invariant under changes of Sq. The NLO impact factor, which is part of 
the scattering amplitude, will change when sq is modified. In our analysis, the NLO impact 
factor was found to be negative and to become more significant when the value of the energy 
scale So is lowered. This is, at least, consistent with the general expectation. 

What remains to be done for the real corrections, is the extension to the case of the transverse 
photon polarisation. It still requires some efforts, as there are additional divergences to be 
dealt with. However, with the tools developed in the longitudinal case and described in this 
paper, we hope to fulfill this task in the near future. 



APPENDIX 



In this appendix we describe our MATHEMATICA and FORTRAN files that contain the 
analytic expressions of our diagrams and are available at: 

http:// www.desy.de/uni-th/smalbc/kyrie/ImpFactor/index.htmll 

A. First, we provide the products of diagrams in momentum space. The >!B's, defined in 
eq.(j2I), are listed in six 2-dim. MATHEMATICA arrays, corresponding to the three groups 
of products, AAij , ABij , BBij and the two color factors Ca,Cf- With AAij = Ca AAfj"^ + 
Cp AA^j' (and correspondingly ABij, BBij, ■ ■ ■) we define the following arrays: 

M_CFaa[[i, j]] = ^^AA^^ , M_CAab[[i, j]] = ^^Sg^ 

and accordingly 

M_CFab[[i,j]], M_CFbb[[i,j]], M_CAaa[[i, j]] , M_CAbb[[i, j]]. 
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The six arrays are stored in the files 

M_CFaa.m, M_CFab.m, M_CFbb.m and M_CAaa.in, M_CAab.in, M_CAbb.m 

In these arrays the following symbols are used: 

L2 = £^ K2 = fc^ R2 = r^ KL = fc£ KR = fcr, LR = £r. 



ep = e, a = a, al = 



Dl . . . D7 are the denominators as given in eqs. (jimifH|) . The BA products are obtained via 
M_CFba[[i, j]] = M_CFab[[j, i]] and M_CAba[[i, j]] = M_CAab[[j , i]]. 

B. We also provide the analytic expressions X, Xgoft, ■ • ■, used in the formulae for the real 
corrections in eqs.lf^ and (j37|) . For any diagram a set of expressions X, Xsoft, • • • is needed. 
In accordance with our labelling of the diagrams, we use subscripts for the X's and for their 
approximations. For instance, [Xsoft]^^-^''"^ corresponds to the Cp part of ABij\soft- We again 
list our results in six MATHEMATICA arrays, but each now having 3 dimensions. The first 
two indices of an array run from 1 to 5 and specifies the diagram. The last index, n, labels 
the approximation; the entry at n=l is the substitution rule for one Feynman parameter as 
imposed by the delta function for that diagram. One is however free to choose any Feynman 
parameter except for the i—uv divergent diagrams where it has to be (see text). Recall 
that in case of these diagrams the substitution of P2 has already been applied to X£„uv If, 
in a given limit, a diagram does not diverge, the corresponding entry in the array is zero. 
The arrays are defined in the following way: 

CFaa[[i, j, n]] = (3^1 — X, Xcen, -^cen,soft, -^soft, -^q, -^q,soft, 

AA,Cf 



^q; ^q,soft; ^i-uv, ~ Arsoft)i 



CAab[[i, j, n]] — P ^ I — Y2P,X, Xcen, -^cen,soft, -^soft, -^q, -^q,soft, 

1 AB, Ca 



, soft 5 



soft ) is 



and, in complete analogy: CFab[[i, j, n]], CFbb[[i, j, n]] and CAaa[[i, j, n]], CAbb[[i, j, n]] . 
The six arrays are available in the files 



CFaa.m, CFab.m, CFbb.m and CAaa.m, CAab.m, CAbb.m 
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The remaining expressions for the BA diagrams are obtained by the replacements: 

CAba[[i, j,n]] = CAab[[j,i,n]] 
CFba[[i,j,n]] = CFab[[j,i,n]]. 

C. In order to calculate numerical values for the real corrections one has to perform the 
integrations over a,a£,l3i which, as indicated in eqs. (jKT|l and (j^Zj), are interrelated by their 
limits of integration. For the numerical computation it is most suitable to have decoupled 
integrals; we therefore introduce new variables Xi such that the integrals extend from to 1 

ae,f3i — >xo,Xi. (A.l) 

For a given diagram the integrand is a sum of expressions, {X — Xgott — ■■■), each of them 
being divergent in some limit. In some cases this sum has to be cast in an analytic form, 
suitable for the numerics. We list here those integrands that we have used for the numerical 
integration. For each diagram we provide the quantities Q, Qj such that the real corrections 
in and (jFTf) read: 



^(l,rea/) 

7* 



finite 



(Svr^ - 28) 



Cp (27r)6 ^24 

+ / da dxQ dxi dx2 dx^ Q'"^ | , (A.2) 



^(l,rea/) 



■7' 



Cf 

Unite e'^ej^N^ - 1 / 7r2 



(37r^ - 2i 



Ca (27r)6 V 48 

+ / da dxQ dxi dx2 dx^ {Q — fti)'"^ | . (A.3) 

Ca / 

The sums extend over the Cp and Ca parts of all diagrams. Q stands for the combination of 
X and its subtractions, as required by eqs. (|51l57|) . expressed in terms of the new variables. 
Qj corresponds to the sixth line in eq.()57|): the substitution of is different from that in Q, 
as a result of the different limits in the ae integration. Note that, in general, the X depend 
on different (3i, and the fi's therefore are functions of different indicated in ()A.1|1 . 

We list the fi's as obtained from the substitution ()A.1|) instead of renaming all variables into 
Xi,X2,X3 as in eqs. ()A.21 IA.3|) . This allows to localize the various limits discussed above. 

The files 



CAintegrand.f 
CFintegrand.f 
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contain, as FORTRAN functions, the f^'s and flj^s for a minimal set of diagrams, separated 
into Ca and Cp parts. All functions have, as arguments, r^, sq, A, a (R2, sO, LA, a), and four 
variables Xi. However, they do not always depend on all of them. The functions are labelled 
in accordance with the diagram they correspond to. In case of Ir-divergent diagrams we 
have combined AAij + ABij {BBij + BAij). 

Labels of the fi's are: 

Aabij : Fabij : AB^^ , 

and correspondingly Aaai j , Abai j , . . . , 

Aaabij : AA^^+AB^^, Abbaij : BB^^+BAf-". 

Labels of the fi/'s are: 

Aaabij I : AA^/" + ABf^^ , , Abbaij I : + i3^g-\ 

The function Faal2(R2, LA, sO, a, xO, xl, x4, x7), for instance, that can be found in the file 
CFintegrand.f has to be inserted in eq. ()A.2|) as in case of the diagram AAi2- AH those 
functions which are not listed in the files CAintegrand.f , CFintegrand.f are either zero or 
can be obtained by exploiting the following symmetries 

Aaaij = Aaaji, Aabij = Abaji, Abbij = Abbji , 

and accordingly for the Cp parts, 
Aaabij + Abbaij = Aaabji + Abbaji, 
Aaabij I + Abbaij I = Aaabjil + Abbajil. 
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